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Abstract: Most recombination-style jet algorithms cluster soft gluons in a complex way. 
This leads to correlations in the soft gluon phase space and introduces logarithmic cor- 
rections to jet cross sections. The leading Abelian clustering logarithms occur at least at 
next-to leading logarithm (NLL) in the exponent of the distribution, and we show that 
new clustering effects contributing at NLL likely arise at each order. Therefore we find 
that it is unlikely that clustering logs can be resummed to NLL. Clustering logarithms 
make the anti-kx algorithm theoretically preferred, for which they are power suppressed. 
They can arise in Abelian and non- Abelian terms, and we calculate the Abelian clustering 
logarithms at 0(u 2 s ) for the jet mass distribution using the Cambridge/ Aachen and kp 
algorithms, including jet radius dependence, which extends previous results. We find that 
previously identified logarithms from clustering effects can be naturally thought of as a 
class of non-global logarithms (NGLs), which have traditionally been tied to non- Abelian 
correlations in soft gluon emission. 
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1. Introduction 

Events with jets are a key signature at collider experiments. Jet cross sections are used to 
measure jet rates and the structure of jets, providing leverage to discriminate between dif- 
ferent signals. Accurate theoretical predictions for jet cross sections require understanding, 
among other things, the perturbative structure of these rates. This is often complicated 
by clustering effects of the jet algorithm, the presence of many scales in the measurement, 
and the need to resum large logarithms. 
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The perturbative series for a jet cross section is governed by two things: the structure 
of QCD matrix elements and the specific observable that is measured. Measurements that 
probe the structure of jets can generate complex perturbative series that depend on many 
kinematic scales. Even basic jet observables depend on the choice of jet algorithm and 
associated parameters, such as the jet radius R and the veto A on soft jets. 

When a measurement is sensitive to soft and collinear radiation, large logarithms can 
appear in the perturbative series and spoil a naive expansion in a s . This occurs for a 
wide range of jet observables of interest. In this paper we focus on the measurement of 
p and A in e + e~ — > 2 jet events, where p is the sum of the two jet invariant masses 
scaled by the center of mass energy Q, p = (m\ + m^/Q 2 , and A is the total energy 
outside of the two jets. The cross section d 2 a/dpdA is sensitive to soft and collinear 
radiation when p Q and has leading logarithms of the form a™ ln 2n p. These logs must 
be resummed to regain perturbative accuracy, and this is achieved through factorization 
and renormalization group evolution (RGE) both in perturbative QCD [1, 2] and in soft- 
collinear effective theory (SCET) [3, 4, 5, 6, 7]. Factorization separates the perturbative 
cross section into separately calculable pieces that depend on different energy scales. For 
the total dijet mass distribution the hierarchy of scales when p <C Q separates the physics 
of the event into distinct processes, given schematically by 



The short distance interaction occurs at a scale of order Q in the hard function H and 
produces the high energy partons that give rise to the jets. This is far above the jet mass 
Q^fp which is the scale associated with the collinear evolution in the jet function J. Soft 
radiation inside of and between jets exists at an even lower energy, Qp, described by the 
soft function, S. 

Non-global logarithms (NGLs) arise when the measurement groups the soft radiation 
in distinct regions of phase space, each of which contributes to a different observable [8, 9]. 
In the total dijet mass cross section, there are NGLs of the scale ratio of the dijet mass 
p and the veto A on additional jets. These start at order a 2 In 2 Qp/ A and continue to 
contribute at next-to-leading log 1 (NLL) at higher orders. NGLs exist due to correlations 
in non-Abelian matrix elements for soft gluon emission. In factorization theorems of the 
form of Eq. (1.1), they enter as terms in the soft function not associated with the anomalous 
dimension. Because they are connected with scales not explicitly separated by the functions 
in the factorization theorem, they are not resummed by renormalization group evolution 



Recombination-style jet algorithms are a common choice to define jets. When perform- 
ing theoretical calculations using these algorithms, the phase space constraints quickly be- 
come very complex as the number of final state partons increases. When a pair of particles 
is clustered, one or both of these partons may be pulled into (or out of) the jet, changing 

1 We will use a standard log counting scheme, counting in the exponent of the distribution. In this case 
N fe LL terms are of order a" \n n+1 ~ k for all n (with k = for leading log). We will also refer to log counting 
in the distribution itself, which sums fewer logs at an equivalent order. 





(RGE). 
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Figure 1: The two particle configurations that change the jet boundary. Particles k\ and k 2 are 
combined across the jet boundary to give the combined pair k t = k\ + k%. Figures (a) and (b) 
show the in-out contribution where particles are clustered across the jet boundary and (c) shows 
the out-out contribution where two particles out of the jet are combined into the jet. 

the jet boundary and producing correlations between particles in the measurement of the 
final state. The schematic configurations for 2 particles are shown in Fig. l(a)-(c). For 
observables sensitive to soft radiation, clustering effects can lead to logarithms in the cross 
section. These are termed clustering logarithms and have been calculated for specific jet 
cross sections [10, 11, 12]. 

In this paper we will show that clustering logarithms arise at least at NLL and that 
they give rise to a new and unique contribution at each order in a s . It is therefore unlikely 
that clustering logarithms can be resummed to NLL, which limits the perturbative accuracy 
of jet observables involving algorithms such as Cambridge/ Aachen (C/A) or kp [13, 14, 
15, 16], for which clustering effects arise. This impacts a wide range of processes involving 
jets both at e + e~ and hadron colliders and makes the anti-kp algorithm ([17]) and event 
shape observables such as ./V-jettiness ([18]) theoretically preferred, for which clustering 
logarithms do not arise. We will also explore the connection between clustering logarithms 
and NGLs in SCET. Both arise in the soft function, and come from correlations between 
soft gluons. The schematic form of the soft function is 



NGLs come from correlations in the squared matrix element, A, which occur only in non- 
Abelian terms. Clustering logs come from correlations in the measurement, A4, and contain 
terms that violate Abelian exponentiation. In Fig. 2, typical configurations that give rise 
to clustering and non-global logs are shown for a dijet event at 0(a^). The two types of 
logarithms share several characteristics: 

• They arise from correlations between soft gluons, in the matrix element for NGLs 
and in the measurement for clustering logs. 

• They are present when the measurement is divided into regions where different phys- 
ical scales constrain the phase space. When correlated soft gluons enter into different 
regions, they produce logs whose argument is the ratio of these scales. 




(1.2) 



- 3 - 




(a) (b) 

Figure 2: Clustering logs and NGLs arise from correlations between soft gluons in separate regions 
of phase space that are sensitive to different physical scales. In (a), clustering logs arise from 
correlations in the measurement due to clustering. In (b), NGLs arise from correlations in the 
non-Abelian matrix element for multiple soft gluon emission. The similarity in origin leads us to 
call both types of logs NGLs. 



• They are associated with soft divergences in the matrix element. 

• They are independent of the UV divergences that give rise to the soft function anoma- 
lous dimension and are therefore not summed by standard RGE. 

We find it natural to consider clustering logs as a class of NGLs, and we will call them 
clustering NGLs. 

We study clustering NGLs in the total dijet mass distribution for e + e~ — > 2 jets 

This observable was defined in [19] and called jet thrust, r^, where uj was used instead of 
A. This distribution depends on the jet algorithm, and in Sec. 2 we give the factorization 
theorem and discuss the impact of different jet algorithms on clustering logs. We will show 
that the anti-kx algorithm obeys Abelian exponentiation and does not generate clustering 
NGLs, but that other iterative jet algorithms do. In Sec. 3 we calculate the O{ol) Abelian 
clustering NGLs for the C/A and bp algorithms away from the small R limit, finding 
excellent agreement with the Monte Carlo program EVENT2 [20, 21]. This extends the 
work of [12], which calculated the leading clustering logarithms at O(a^) in the small R 
limit for the C/A algorithm. 

In Sec. 4 we discuss the properties of clustering NGLs and their connection to tradi- 
tional NGLs, and in Sec. 5 we discuss the structure of Abelian clustering NGLs at higher 
orders. Since the Abelian matrix element is simple at any order, we need only understand 
how the measurement function for a jet algorithm handles an n particle final state. We 
show that at (D(a™) Abelian clustering NGLs contribute at NLL order, a™ ln n (Qp/A). Since 
this is the first order at which n particle clustering effects arise, we expect there will be a 
new contribution to the coefficient of NLL clustering logarithms at each order in a s that 



- 4 - 



is unrelated to lower order coefficients. We see this explicitly from our results at 0(a 2 ). 
Therefore, it seems that clustering NGLs cannot be resummed, since determining the cross 
section at NLL would require the calculation of an infinite number of distinct coefficients. 
In [22] we prove that this is in fact the case by using a novel framework to express the 
all-orders form of the measurement function. However, NLL resummation using weaker 
log counting in the distribution itself is still valid, as this counting reorganizes the terms 
to make the lowest orders (in a s ) of clustering NGLs the most important. In Sec. 6 we 
present our conclusions. 



2. The Dijet Mass Distribution 

The process e + e~ — > 2 jets provides a simple system to study the perturbative structure of 
jet cross sections. We focus on the total dijet mass distribution Eq. (1.3), which is obtained 
as follows. For each event, we cluster the final state into jets using a jet algorithm and label 
the invariant mass of the two most energetic jets mi and mi. We measure p = (mf+m 2 ,) /Q 2 
and the total energy A of particles outside the two leading jets. Additional jets are vetoed 
by requiring A <^ Q. 

Adding the jets' masses together reduces the complexity of the calculations by removing 
one scale from the problem. The doubly differential cross section, Eq. (1.3), and related 
observables have been used to study various properties of jet cross sections, including the 
cross section dependence on the jet algorithm, the jet radius R, the ability to simultaneously 
resum logarithms of ratios of the hard scale Q, the jet scale Qy/p, and the soft scales Qp 
and A, and the impact of non-global logarithms [23, 24, 10, 11, 25, 26, 27, 28, 29, 30, 19, 31]. 

We will study this cross section using the framework of SCET, as well as the double 
cumulant when convenient, 

f pc f Ac d 2 a 
Efc,A ) = / o i P j o dA — . (2.1) 

We will work in the dijet limit of small jet mass and small out-of-jet energy, 

{Qp,A}«Q^«Q, (2.2) 

where Q ~ Ej and the final state is dominated by soft and collinear radiation. SCET can 
be applied in this regime, and the cross section factorizes into hard, jet, and soft functions 
[32, 33, 34]: 



d 2 a 
dpdA 



a H 2 (Q) J dp n dp n dp s 5(p - p n - p n - p s ) Jn(Pn) Jfi(Pn) S(p s ,A) , (2.3) 



where the labels n and n refer to the back-to-back jet directions 2 ra^ra^ = (1, ±n). Each 
of these functions describe the physics at different scales, depicted in Fig. 3. The hard 



2 Any four-vector can be decomposed as 



u " _ W u 

P = n • p-j- + n ■ p— + pi 
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Figure 3: Relevant scales in the cross-section with a measurement on the sum of jet masses with 
a veto on the radiation outside the jet. The soft function depends on both the scales p™ ~ Q p and 



p° s ut ~ A. 



function describes the short distance interaction at the hard scale Q and is independent of 
the observable. The jet function describes collinear evolution of massless partons into jets 
with mass m ~ Q J~p, and the soft function describes global soft radiation in and between 
jets. The radiation described by the jet function is collimated along the jet axis. As long as 
this collinear radiation is sufficiently narrow compared to the jet radius, satisfying p -C R 2 , 
it will be insensitive to the jet boundary and therefore different jet algorithms will cluster 
these particles equivalently [25, 26, 35]. In this limit the jet functions J n ^ are universal 
up to power corrections for a wide range of processes. Soft radiation, however, probes the 
entire final state, and is sensitive to scales in the jet (the scale Qp) and out of the jet (the 
scale A, as well as the parameter R). It is therefore sensitive to the jet boundary and the 
soft phase space constraints will differ between jet algorithms, leading us to focus on the 
soft function. 

The soft function is a. forward, scattering matrix clement of soft Wilson lines Y n Y^ 

S(p, A) = ~TY <0| Y n y n t M(p, A) Y n Y*\0) , (2.4) 

where the measurement operator M (p, A) implements the measurement on the final state, 
including the jet algorithm. The measurement acts on the final state of soft particles with 
momenta {k s } as follows 



M(p, A) \X S ) = M (p, A; {k s }) \X a ) , (2.5) 



where 



M{p,A;{k s })=6 (p-i n -^-\> E ) 6 A " E k i) ■ ( 2 - 6 ) 

iSn jet idn jet / \ i £n,n jets 
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The jet algorithm determines which particles are in the jets, meaning the sums over particles 
in and out of the jets have phase space constraints imposed by the jet algorithm. 

Soft radiation will affect the value of the jet mass, and thus the observable p; hence, 
the soft function is convolved with the jet function in the factorization theorem of Eq. (2.3). 
It will be convenient to express the cross section in Fourier space for both p and A, which 
removes the convolutions and makes the measurement more straightforward: 

d 2 a f°° d 2 a 

/ dpdk- — — - exp(— ixp) exp(— iyA/Q) 
y_oo dpdA 

a H 2 (Q)J n (x)J n (x)S(x,y). (2.7) 



dx dy 



The Fourier space jet and soft functions, J and S, are the Fourier transforms of the 
momentum space functions. The Fourier space form of the measurement in the soft sector 
is 

M(x,y;{k s }) = exp ( -ix ^ ^ — J exp | -ix ^ ^— — J exp | -iy ^ ^ 

\ iSn jet / \ iGn jet / \ i ^n,n jets 

(2.8) 

This form is useful as it removes the convolution in Eq. (2.6) between different particles' 
momenta, and will be helpful in understanding NGLs. 

2.1 Jet Algorithms and Boundary Clustering 

Jet algorithms are used to identify energetic clusters of radiation. Inclusive recombination 
algorithms, which we study here, build the jet through pairwise clustering of particles. 
These algorithms utilize two kinds of metrics; a pairwise metric pij between two particles 
and a single particle metric pi for an individual particle. The algorithm works recursively, 
finding the minimum of all pairwise and single particle metrics. If the minimum is a 
pairwise metric, then that pair is merged by adding their four momenta, "Pi+j = Pi +Pj- 
If the minimum is a single particle metric, then that particle is promoted to a candidate 
jet. This is repeated until all particles have been clustered into candidate jets. Finally, an 
energy or pt cut is used to veto soft jets and select final state jets. 

The kx class is a common set of jet algorithms [13, 14, 15, 16, 17]. The metric for 
these algorithms is parameterized by a number a. For e + e~ collisions, the metrics take the 
form: 

Pij = 2 min {Ef a , Ef) (1 - cos By) , 

Pi = 2Ef a {l - cos R) , (2.9) 

where Oij is the angle between particles i and j. The radius R sets the size of the jet, and 
is the maximum angle between two particles for which a single clustering can occur. This 
class includes the kx, Cambridge/ Aachen (C/A), and anti-kp algorithms: 

a = 1 : k T , a = : C/A , a = -1 : anti-k T . (2.10) 
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For a single soft particle, these algorithms give the same phase space constraints; the 
jet boundary is simply a cone of radius 22 around the jet direction. The position space 
measurement function for a single particle for kp, C/A and anti-kr is 

M^\x,y-k) = exp(-ix^@ k n -ix^e k n -iy^@ k on ^ , (2.11) 

where 

k 6 n jet : 9* = 0(22 - 6 kn ) , 
k e n iet : Q k n = 6(9 kn - (tt - R)) , 

k$n,fi jets : 9* ut = 0((tt - R) - 6 kn ) 6(e kn - R) , (2.12) 

In the small p and small A regime, when p <C R 2 , collinear radiation in the jet is 
parametrically narrower than the jet size and is insensitive to clustering effects of the 
algorithm. The jet boundary is therefore determined by soft radiation. The behavior of 
the kx class of algorithms have been previously analyzed in the context of SCET, where 
the metrics determine the characteristic sequence of clustering and the effect on the jet 
boundary [36]. Soft partons will either first cluster with collinear radiation or among 
themselves, depending on the algorithm. In the former case, the phase space restrictions 
take the form of a cone of radius R around the jet axis. When soft particles merge among 
themselves, they can change the geometry of the jet boundary away from a cone. We refer 
to this as 

boundary clustering : clustering of soft particles across the jet boundary of radius R 
around the jet direction. 

While clustering purely inside or outside of the jet of course still occurs, it does not change 
the measurement function, Eq. (2.8), and therefore we do not need to consider it. Boundary 
clustering however is relevant because it leads to a change in the observable. As an example 
consider the configuration with two nearby soft gluons shown in Fig. 1(a) and (b), where 
gluon 1 is inside the cone of the n jet {9\ n < R) and gluon 2 is outside this cone (#2n > 22). 
If these gluons are not clustered, then gluon 1 is in the jet and contributes to the p 
measurement and gluon 2 is out of the jet and contributes to the A measurement. However, 
if the gluons are clustered, then both gluons will either be in the jet, contributing to p, or 
be out of the jet, contributing to A, as shown in Fig. 1(a) and Fig. 1(b), respectively. Thus, 
boundary clustering changes the observable and can introduce logarithmic corrections in 
the soft regime. 

Boundary clustering gives correlations in the phase space constraints for soft particles. 
It is useful to express the measurement function in Eq. (2.8) for n particles in terms of a 
product of independent single particle constraints plus a correction A.M a i g due to boundary 
clustering: 

n 

M%l(x, y, {h,..., k n }) = Y[M^(x, y, h) + AM alg (x, y,{h,---, K}) . (2.13) 

8=1 
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Figure 4: The region of soft particle boundary clustering is shown schematically in grey for 
(a) anti-kx and (b) C/A, kx- This corresponds to the correction to a cone of radius R in the 
measurement function given by AA4 a i g for each algorithm. For anti-kx this region of phase space 
is power suppressed by a factor A 2 ~ p <C 1. For C/A and kx the region of soft particle boundary 
clustering corresponds to an 0(1) region of phase space. 

In Fig. 4, we show the schematic size of the region contributing to A.A/f a i g for the anti- 
kx, C/A, and kx algorithms. Since the anti-kx metric is weighted by the inverse energy, 
in the dijet limit the algorithm first clusters particles with scaling collinear to the jet 
directions. Then clustering between soft and collinear radiation typically occurs before 
clustering among soft particles, meaning that the jet boundary for soft particles is only 
set by the jet axis and is a cone of radius R up to power corrections [36]. As a result, for 
the anti-kx algorithm, the region of phase space where boundary clustering can take place 
(AA4 a k T ) is power suppressed and the measurement function in Eq. (2.13) factorizes: 



The measurement function also factorizes in this way for global event shape observables 
such as thrust, angularities, or A^-jettiness [32, 34, 37, 18]. For the C/A and kx algorithms, 
boundary clustering can change the jet boundary by an 0(1) amount and therefore con- 
tributes at leading power. Since the C/A metric depends only on angle, pairs of soft par- 
ticles that are closer to each other than the jet axis will cluster. For the kx algorithm, this 
clustering is enhanced, since the metric is weighted to preferentially merge softer particles 
earlier in the algorithm. 

The correction term AA4 a i g in the measurement function leads to a correction in the 
soft function, 



n 




(2.14) 



i=i 



S & \ g (x, y) = S akT (x,y) + AS a i g (x, y) 



(2.15) 
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and correspondingly in the cross section, 



dx dy dx dy dx dy 



(2.16) 



Boundary clustering will generate clustering logs in the cross section. We will focus on the 
Abelian terms, where the analysis is simplest. 

2.2 Abelian Exponentiation 

At C(a"), a general Fourier space soft function contribution has the form: 



S<">({*}) 



n 



d% 



A^[{k s })M^[{x},{k s }) 



(2.17) 



where A^ is the n-loop squared matrix element and implements the measurement 

of the Fourier space observables {x} on the n-particle final state. Abelian matrix elements 
for n particles factorize: 



1 n 

A (n) {{k s }) = ^7 II" 4 (^)> 
where the one loop matrix element is 



(2.18) 



8=1 



1 



n ■ ki)(n ■ ki 



(2.19) 



For observables for which the measurement function factorizes, as in Eq. (2.14), the Abelian 
contributions to the soft function also factorize. This implies Abelian exponentiation, where 
we use the common definition that the n-loop Abelian contribution to the soft function is 
determined solely by the one-loop contribution: 



sSel.(M) 



1 

111 
1 

n! 



d A k 



(2vr) 4 
£«({*}) 



A w (k)M {1) ({x},k) 



which implies 



SAbci.(M) =exp \sW({x}) 



(2.20) 



(2.21) 



The renormalized one-loop soft function in Fourier space is (for all algorithms in the 
bp class) 



S W (x,y) 



o. s Cf 1 
vr Q 



f ^ tan | ie 7£ xj + 2 In (ten 2 | ) — 



R 



2Q 



le ' y 



- In 2 (tan 2 — 
2 V 2 



2Li 2 



tan 



R 



57T 2 

~12 



(2.22) 
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Figure 5: A configuration of two soft gluons that will boundary cluster. Even though gluon 1 is 
close to the jet axis, gluon 2 is sufficiently soft so that it will cluster with gluon 1 before gluon 1 
clusters with the jet. 

2.3 Violation of Abelian Exponentiation 

Corrections in the soft function due to boundary clustering, AS'aig, violate Abelian exponen- 
tiation, since the n-loop Abelian contribution to the soft function is no longer determined 
by the one-loop contribution. This effect starts at O(a^) and the relevant contributions 
at this order are shown in Fig. 1. Such configurations give logarithms that are sensitive 
to both p and A in momentum space, as we shall see by explicit calculation at O(a^Cp) 
in Sec. 3. In general each final state gluon can contribute a double logarithm to the cross 
section, a" ln 2n (Qp/A), with the double log associated with soft and collinear divergences. 
Unlike the non-Abelian case, the 0(a™) Abelian soft matrix element, Eqs. (2.18) and (2.19), 
has no l/k{ ■ kj structures and so the only collinear divergences in the matrix element arise 
for soft gluons collinear to the jet directions. 

However, in order for gluons to boundary cluster and contribute to AS a i g they cannot 
be collinear to the jet direction. For C/A the closest a soft gluon can be to the jet axis 
and be merged across the jet boundary is 2~( n-1 )i? for n final state gluons. For k/r, the 
condition is more subtle. Suppose we have two soft gluons 1 and 2, with gluon 1 at a small 
angle 9\ to the jet axis, and gluon 2 at a larger angle 9% to the jet axis (9± <C O2). If the 
two soft gluons cluster instead of gluon 1 clustering with the jet, then 

9 9 

E 2 sin -y < Ei sin -± . (2.23) 

Such a configuration is shown in Fig. 5. Since gluon 1 is close to the jet axis, 9\ <C #12, 
meaning gluon 2 must be much softer than gluon 1: 

^1 < ^i/ 2 (2 24 ) 

£1 sin 012/2" 1 > 

This means that the combined pair will be nearly collinear with gluon 1. If gluon 2 is 
outside of the jet, such that the clustering contributes to AS'aig, then since #12 < R we 
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must have 

R<6 2 <R + 6i. (2.25) 

This is a narrow range of angles, meaning the phase space for gluon 2 is power suppressed. 
Although these configurations will contribute to the clustering logs at leading power, the 
power suppression of phase space for gluon 2 mitigates the collinear logarithmic enhance- 
ments from soft gluon clustering with kp. Since the momentum of a soft gluon close to the 
jet axis is nearly unchanged if it clusters with another soft gluon instead of the jet, this 
argument extends to configurations with more soft gluons. 

For the Abelian terms for both the C/A and kp algorithms, we see that if a soft gluon 
contributes to AS a \g(p, A) then it can contribute at most a single log from a soft divergence, 
since collinear divergences with the jet direction are excluded by boundary clustering. This 
provides a very straightforward log counting of the Abelian terms in AS'aig: for n final state 
gluons, the leading term in AS^ is at least of order a" ln n (Qp/A), which is NLL. 

3. Abelian Clustering Non-Global Logarithms for the Dijet Mass Distri- 
bution 

In this section we calculate the 0(a 2 ) violations to Abelian exponentiation from soft clus- 
tering. We focus on the C/A and kp algorithms, but the same effect occurs for any jet 
algorithm that does not satisfy Eq. (2.14). It is useful to express the measurement for a 
generic algorithm in terms of the anti-kp measurement function .M a k T and a correction 
term as in Eq. (2.13). This gives a corresponding correction to the soft function AS^g, as 
in Eq. (2.15), which at 0(a 2 ) can be written as 



q( 2 ) - 1 

^alg " 2 



D ak T 



+ A^2- (3.1) 



The effect of correlations in the soft gluon phase space are expressed entirely by the soft 
function correction term. In this section we work in momentum space, where analogous 

(2) 

expressions hold, and we calculate AS^(p,A) for C/A and kp. 
3.1 Abelian Non-Global Logarithms at 0{a 2 s ) 

The Abelian soft function correction term at 0(oq) due to clustering with an algorithm is 
ASig(p,A) = \ J ^^A^Hk^K^AM^MhM), (3-2) 

where A^ l \ki^) is given in Eq. (2.19). Note that in pure dimensional regularization, a 2 C\ 
terms with a virtual gluon are scaleless and hence zero - these terms make the soft function 
IR finite. 

The configurations which contribute to the correction terms in the measurement func- 

(2) 

tion, AMgjg, are shown in Fig. 1. The solid line represents the anit-kx jet boundary of 
radius R around the jet direction. The clustering of soft gluons in and out of the boundary 
is represented by Fig. 1(a) and (b), which we refer to as the in-out contribution (10), while 
Fig. 1(c) shows two gluons out of the jet being clustered in to the jet, which we refer to 
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as the out-out contribution (OO). It is useful to define the constraints requiring the soft 
gluons to be in a region of phase space where the algorithm produces a different cross sec- 
tion than the anti-kx algorithm as @f g and ©2j? for the in-out and out-out contributions 
respectively. We can write these schematically as 

®aig = 0(^1 m n J e t) 9(k 2 out of jets) 9(k\, k 2 cluster) 

8 al g = 6{k\ out of jets) 6(k 2 out of jets) 9(k\, k 2 cluster) (3.3) 

For each contribution Fig. l(a)-(c), we must remove the contribution of the anti-kx mea- 
surement to that region of phase space included in A^ a k x in the momentum space version 
of Eq. (2.13). We can write the measurement function correction term as 

AM ( V(p, A; kl , k 2 ) = 4 eg [Mf + Mf + Mf) + 2 eg° [M° p ° + M° a °] , (3.4) 

where the in-out contributions from Fig. 1(a) and (b) are 

Mf = 0(h + k 2 mn jet) S(p - (k+ + k+)/Q) 5(A) , 

Mf = 0{ki + k 2 out of jets) 6(p) 5(A - (fc? + k° 2 )) , 

Mf = -S(p- kf/Q) 5(A - k° 2 ) . (3.5) 

Mf is the contribution where the algorithm clusters the gluons into the jet, Mf is the 
contribution where the algorithm clusters the gluons out of the jet, and Mf removes the 
contribution from the anti-kx algorithm in the correction term. The out-out contributions 
from Fig. 1(c) are 

M° p ° = 9(ki + k 2 in n jet) S(p - (kf + k+)/Q) 5{A) , 

M°° = - 9(ki + k 2 mn jet) S(p - kf/Q) 5(A - k%) . (3.6) 

As with the in-out contributions, M® is the contribution where the algorithm clusters 
the gluons into the jet and M®° removes the contribution from the anti-kx algorithm. 
Note the combinatoric factors in the measurement function arise because interchanging k\ 
and k 2 produces an identical contribution to the cross section and since there are two jets, 
ki can be in either jet for the in-out contribution. 

We present the in-out regions of phase space for the C/A algorithm schematically in 
Fig. 6. The region in red corresponds to gluon 1 in the jet, where boundary clustering can 
take place, while the region in blue contributes to the measurement of the jet mass cumulant 
p c . The overlap between these two regions corresponds to Mf, the region outside up to 
k± + /c!> < A contributes to Mf and the combined regions are subtracted in Mf. For 
the Abelian terms, collinear divergences only arise along the axes, while soft divergences 
occur at the origin. It is clear from the figure that boundary clustering (red region) does 
not contain any collinear divergences and that the while each term in Eq. (3.5) contributes 
soft divergences these are removed in the sum. We can see this explicitly by taking the 
soft limit of M^ g in Eq. (3.4), which vanishes as both k^ — > and k® — > 0. Fig. 6 also 
illustrates that the jet mass measurement (blue region) constrains the phase space in the 
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C/A phase space for in-out boundary clustering 



k 



+ t 




< pcQ 



k: 



Figure 6: The C/A in-out region of phase space at 0(a 2 s ) is shown schematically, where we overlay 
the phase space for gluon 1 and 2 in the (k^~ , k~) plane. The region where the gluon in the jet can 
boundary cluster is shaded in red (dark gray). The region contributing to the jet mass cumulant 
p c is shaded in blue (light gray). The over lap between these two regions corresponds to A4 l ®, 
the region outside the jets up to k\ + k® < A corresponds to A^Jf and the combined regions 
are subtracted in A4 l ®, removing the IR divergences from this contribution to the measurement 
function. 



UV rather than the jet boundary R and therefore the soft function anomalous dimension 
is independent of the algorithm. 

We present the results of the calculations of Eq. (3.2) for C/A and k/p here and leave 
the details to Appendix A. As we have argued, the logarithms in AS , a i g (/j, A) come from 
soft divergences, which we extract by integrating over and k®. We express the final 
result in terms of coefficients of the double NGL (C^(R), which comes from in-out terms 



only), and the single NGL (C$^(R) + C§£ 



alg 

7jq a ig(-R) + Coo alg(^))> wn i° n we calculate numerically and 
plot in Fig. 7. The integrals that determine these coefficients are given in Eqs. (A. 45) and 
(B.15). For the bp algorithm, these coefficients depend on the values of p and A. This 
makes it difficult to write AS^ in terms of standard distribution functions, and so for the 
kp algorithm we only determine the cumulant. 

For the C/A algorithm, the correction terms to the soft function due to correlations 
introduced by clustering are 
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Abelian clustering NGL coefficients for C/A 
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Abelian clustering NGL coefficients for kj 
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Figure 7: The coefficients of the in-out double NGL, C^(R), the in-out single NGL, C$ alg (R), 

and the out-out single NGL, Cqq alg (R), as a function of R, for C/A (left) and kx (right). These 

arc the coefficients of (a s CV/7r) 2 In 2,1 Qp/2Asui 2 § in the cumulant AX3g . For kx, the in-out 
double and single log coefficients depend on the value of the argument of the NGLs, and we plot 
the coefficients for Qp/2Asin 2 ~ > 1. They are constant over this range. 



+ (constant) 5 
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2 sin 
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Cq and L\ are standard plus distributions [38], 

~0{x) 



£ (x) 



(x) Inx 



(3.7) 



(3.8) 



The term proportional to 5(p)5(A) is a number dependent only on R that we do not 
calculate. For the cumulant we find the same form for both algorithms, 



AS 



(2) 



Qpc 



2A C sin 2 f 



+ constant 



(3.9) 



For the bp algorithm the in-out single and double log coefficients depend on the value of 
the argument of the log, Qp c /2A C sin 2 R/2. As discussed in Appendix A, when this ratio 
is larger than 1 the coefficients are constant, and decrease as this ratio becomes smaller 
than 1. In Fig. 7, we plot the kx coefficients for the regime Qp c /2A c sm 2 R/2 > 1. In the 
next section we compare our results to the distribution calculated using the Monte Carlo 
program EVENT2. 

The coefficient of the leading double logarithm of o? s C\ In 2 p c in the correction term 
for the cumulant was calculated for C/A in the limit of small jet size, R, in [12]. We find 
our result for the double log coefficient for C/A in the limit of small R, 



,(0) L (i?) = 0.183, 

agrees within the precision of the value in [12], which is 0.185 in our units. 



(3.10) 
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3.2 Comparison to EVENT2 

We can compare the results of our calculations to the output of the Monte Carlo program 
EVENT2 [20, 21]. EVENT2 implements Catani-Seymour subtractions and contains the 
matrix elements necessary to compute observables in e + e~ events to 0{oq) that vanish 
in the 2-jet limit. We use EVENT2 to calculate the dijet mass distribution, taking the 
cumulant in the out-of-jet energy cut A c : 



EVENT2 calculates a binned distribution in p, with the cross section in a bin [p m i n , Pmax] 



where S is the double cumulant of the cross section, Eq. (2.1). 

We compare the O(a^Cp) terms to EVENT2, first subtracting the global logs from the 
distribution from EVENT2. These global logs, equal to the anti-kx cross section, are given 
by the universal hard and jet functions and the anti-kx soft function through 0(a 2 ); for 
reference, the global logs are given in Section V of [30]. We compare the analytic results 
for the anti-kx, C/A, and kx algorithms to EVENT2. The anti-kx comparison serves 
to verify that the global logs are correct, and the C/A and kx comparisons allow us to 
directly test the calculations of the Abelian NGLs. For each algorithm, we use EVENT2 
to calculate 0"ev2 for A c /Q = 10 -5 and R = 0.2, 0.6, and 1.0. We generated approximately 
2 • 10 12 events for the anti-kx and C/A algorithms and approximately 8 • 10 12 events for 
the kx algorithm. The kx results are less numerically stable, but still allow us to test the 
calculation of the NGLs. 

In Fig. 8, we show the distributions from EVENT2 that confirm the calculations of 
the clustering NGLs. For the anti-kx algorithm, we show the distribution da /dp before 
and after the global logs have been subtracted. After the global logs are removed, the 
distribution vanishes for small p, indicating that only power suppressed terms in p remain 
and verifying the global logs. For the C/A and kx algorithms, we subtract the global 
logs and show the distribution da /dp before and after the clustering NGLs have been 
subtracted. Before the clustering NGLs are subtracted, the distributions for small p have a 
non-zero slope, indicating double logs remain. After they are subtracted, the distribution 
vanishes for small p, confirming our calculations. 

4. Non-Global and Clustering Logarithms 

In this section we discuss general properties of clustering logs that we saw arise explicitly 
in the C(a 2 ) calculation in the previous section. Clustering logs share many of the same 
properties as traditional NGLs, leading us to think of them as a class of NGLs (clustering 
NGLs). We will discuss the points of similarity outlined in the introduction in further 
detail. 




(3.11) 



given by 




dp—(A c ) = S(p max , A c ) - S(/) min ,A c ) 



(3.12) 
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k T algorithm 




Figure 8: Numerical analysis of the EVENT2 results with A c /Q = 10~ 5 . In (a), the Abelian terms 
in the distribution da /dp for the anti-kx algorithm are shown. The distribution after the global 
logs have been subtracted is shown in (b). The result tends to for small p, verifying the global 
logs. In (c) and (e), the distributions for the C/A and kx algorithms after the global logs have 
been subtracted are shown. The distributions are not flat for small p, meaning double logarithms 
remain. These are the clustering NGLs, which are subtracted in (d) and (f). The distribution goes 
to for small p, confirming the calculations of the clustering NGLs. 
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In order for clustering logs to arise, two properties must be satisfied. First, there must 
be correlations in the phase space constraints on multiple gluons in the final state, requiring 
that A.M in Eq. (2.13) is non-zero. Second, the measurement must be divided into regions 
with different scales contributing to the cross section, namely AAA must depend on multiple 
scales. These conditions are satisfied by boundary clustering, which requires at least one 
particle to move between two measurement regions and gives rise to clustering logs of the 
ratio of scales associated with these regions. 

The conditions required for traditional NGLs to arise are completely analogous. There 
must be a correlation between the soft gluons from the non-Abelian matrix elements and 
there must exist a boundary in phase space where the measurements differ on opposite 
side of the boundary. The latter is a feature of the observable, and many jet observables of 
interest contain NGLs. A typical configuration that gives rise to NGLs is shown in Fig. 2(b). 
NGLs have been studied for a variety of observables at 0(a 2 s ). For the distribution we study 
here, d 2 a/dpdA, the non-Abelian NGLs are 3 



The correction due to clustering in the measurement function, AA^aig; gives rise to a non- 
Abelian clustering NGL with a coefficient A/ a i g that depends on the jet algorithm. For 
the anti-kx and cone algorithms, for which clustering logarithms do not contribute, A/ a k T 
is power suppressed. The coefficient / a k x can be determined analytically [30, 19], and is 



The contribution of non-Abelian clustering NGLs, A/ a i g , has been calculated for C/A and 
bp numerically in [30]. These reduce the magnitude of the coefficient in Eq. (4.1) for 
algorithms where soft gluon clustering plays a role. 

We saw from the the argument in Sec. 2.3 and the explicit calculation in Sec. 3 that 
clustering NGLs arise from soft singularities. This is true for both Abelian and non-Abelian 
clustering NGLs. When boundary clustering merges a particle in the jet with one out of 
the jet, the particle in the jet cannot be collinear to the jet axis and hence cannot be 
sensitive to the collinear singularity. This is required by a collinear safe jet algorithm: if a 
soft particle is collinear with the jet axis, it must cluster with the jet. At leading power, 
there will be some region around the jet axis where boundary clustering cannot take place. 
This implies that clustering NGLs arise from soft divergences rather than collinear. 

Traditional NGLs also arise from soft divergences. For non-Abelian real emission terms, 
in addition to the singularities corresponding to the partons becoming soft or collinear to 
the jet direction, there is an additional collinear singularity when the two partons align 
with each other. This singularity generates logarithms of a single scale if the the partons 
propagate into the same region, or two scales if they propagate into different regions. 

3 The power of tan 4r in Eq. (4.1) is different from that in [30]. Since we are not working in the small 
R limit, both represent valid choices. With the scaling of R included in the argument of the NGLs in this 
paper, the coefficients of the a 2 In 2 and a 2 In terms are finite as R — > 0. 




(4.1) 




(4.2) 
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The relationship between these terms contributing to the NGL was explored in [39, 30]. 
For configurations with partons in different measurement regions (such as Fig. 2(b)), the 
collinear singularity is never realized due to the boundary. A collinear enhancement remains 
in the matrix element, and the region of phase space with the two soft partons near the 
boundary dominates the value of the leading NGL coefficient. 

Finally, we can show that clustering NGLs are not associated with any UV or IR 
divergences in the total soft function. The lack of IR divergences follows directly from 
the fact that anti-kx, C/A, and kp are IR safe jet algorithms. The difference in their 
contribution to the soft function, and the cross section, cannot contain any IR divergences. 
Since boundary clustering does not occur for the anti-kp algorithm, boundary clustering 
cannot give IR divergences for any IR safe algorithm. The lack of UV divergences follows 
from the consistency condition on the factorization theorem in Eq. (2.3). Consistency 
requires 

ih{Q) + jjAVpQ) + Un(VpQ) + is(pQ) = o . (4.3) 

For the dijet mass distribution, when p <C R 2 the jet function is inclusive and is indepen- 
dent of the choice of algorithm. This is true for other observables where the phase space 
constraints in the far ultraviolet are more restrictive than the jet algorithm [40]. In this 
case, the hard function and jet function anomalous dimensions are independent of the jet 
algorithm, and Eq. (4.3) requires that the soft function anomalous dimension is indepen- 
dent of the jet algorithm as well. As far as RGE is concerned, there is no difference between 
the anti-kp, C/A, and kp algorithms. This means that clustering effects cannot contribute 
to the soft anomalous dimension, and hence cannot contribute UV divergences to the soft 
function or the cross section. 

Similarly, traditional NGLs are independent of the UV divergences in the soft anoma- 
lous dimension. Since the hard and jet function anomalous dimensions are independent of 
the scale A, Eq. (4.3) implies that the soft anomalous dimension is independent of A as well. 
The non-Abelian NGL involving the scale A in Eq. (4.1) must therefore be independent 
of the renormalization scale p and not contribute the soft function anomalous dimension. 
The details of this structure are described in [30]. 

5. Non-Global Logarithms at All Orders 

Having confirmed the existence of clustering NGLs and established their general properties 
and connection to traditional NGLs, we discuss the general structure of all NGLs beyond 
O(o,g), focusing on the implications for resummation. Resummation enforces a relationship 
among the coefficients in the perturbative expansion of the soft function at a given order 
in an a™ ln m counting scheme. Consider for example the perturbative structure of the 
cumulant soft function. The perturbative series can be arranged according to log counting 
in the exponent as follows 

T, s = exp [lg (a s l) + gi(a s L) + a s g 2 (a s L) H ] (5.1) 

where L represent logs of p or A. The LL series is given by L go(a s L) = d\2 a s L 2 + <i23 a^L 3 -!- 
^34 a 3 L 4 + • • • , the NLL series is given by gi(a s L) = di±a s L + d 2 2 a 2 V. 2 + ■ ■ ■ and the NNLL 
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series by a s g2(a s \-) = ^io a s + <^2i a 1 L + • • • • Resummation at order N fc LL in the exponent 
allows all gi for i < k terms to be included. 

Since clustering NGLs arise due to correlations in soft gluon phase space introduced 
by the jet algorithm, the specific structure of clustering logs at higher orders depends on 
the algorithm's metric and the number of gluons in the final state. As we have seen, the 
kx class of jet algorithms have the same single soft gluon phase space constraints and 
therefore agree at 0{a s ). However, they disagree when two or more gluons are in the final 
state, starting with the coefficient d,22 in the NLL series for both the Abelian and non- 
Abelian terms. We expect this disagreement to continue at higher orders. For example, 
the clustering constraints at O(af) include the region 

' — (3) 

A.Mgjg D (gluons 1,2 and 3 boundary cluster) , (5.2) 

which arises first at this order and is unrelated to the clustering effects at 0(a 2 ) for 
two soft gluons. As demonstrated in Sec. 2.3, the constraint of clustering means this 
contribution will be 0(afL 3 ), thus it adds to the ^33 term in the Abelian NLL series. 
Boundary clustering of n gluons prevents collinear divergences with the jet direction; thus 
Abelian clustering NGLs arise at order a™L n . When only a subset of the soft gluons are 

clustered, the Abelian clustering NGLs could enter at order a™L m , where m > n. As a 

— (3) 

good example, again consider the AjW^ correction, which also includes the region 

' — (3) 

D 9 (gluon 1 is in the jet and does not boundary cluster) 

x 9 (gluon 2 and 3 boundary cluster) . (5-3) 

This will contribute at order aj?L 4 to the Abelian terms because gluon 1 is sensitive to both 
soft and collinear divergences and so can contribute a s L 2 to the cumulant, while gluons 
2 and 3 are sensitive to soft divergences only and therefore contribute a 2 L 2 . This will 
provide a new contribution from clustering to the LL series coefficient (I34 unless the phase 
space constraints in Eq. (5.3) can be expressed in terms of lower order constraints. That 
is, if gluon 1 gives di2a s L 2 and gluons 2 and 3 give the same contribution from boundary 
clustering at 2 loops, c^c^L 2 , the clustering effects from Eq. (5.3) will still be NLL. This 
is non-trivial to show and we leave the proof to future work [22], where we use a novel 
framework to express the all-orders form of the measurement function and use this to show 
that Abelian clustering NGLs contribute only at order a"L n for all n > 2. 

This presents a troubling picture for resummation of clustering NGLs. Unless the 
phase space constraints giving the leading clustering NGLs for n > 3 particles are related 
to lower order phase space constraints for some algorithm, the coefficients in the NLL 
series in the exponent will be unrelated. There is no a priori reason to believe that such a 
property exists, and considering the new Abelian clustering contributions at a 2 and from 
Eq. (5.2) at af, such a relationship between phase space constraints seem very unlikely. 
This is further supported by the general structure of the measurement function that we 
prove in [22]. We therefore find that the NLL resummation of clustering NGLs is likely 
impossible. This is not necessarily true for traditional non- Abelian NGLs, as we discuss 
shortly. 
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We note that, even though resummation is likely broken at NLL, we might be interested 
in log counting in the distribution instead instead of the exponent. In this case, NLL means 
resumming all terms through order a™ln 2ri_2 ; therefore, it is sufficient to perform leading 
log resummation in the exponent and then add the fixed order contribution from the 
clustering NGLs. 

Non-Abelian NGLs arise from correlations in the matrix element for multiple soft gluon 
production. These NGLs are determined by the perturbative structure of matrix elements 
of the soft function, which is the eikonal limit of QCD. Unlike NGLs from soft gluon clus- 
tering, there is reason to believe that NGLs at all orders may be related. The fundamental 
problem with resumming clustering NGLs is that the action of the jet algorithm at higher 
orders involves phase space constraints that are not present at lower orders. Eikonal matrix 
elements at higher orders in QCD are related to lower order matrix elements, especially 
when organized by log counting in the exponent of the soft function. 

Consider an observable, such as thrust (r), that does not have NGLs. Like the anti-kp 
algorithm, the measurement function for thrust is a product of single particle constraints. 
Resummation can be performed in the exponent of the distribution in this case. For 
thrust, this has been performed to N 3 LL [41]. For observables of this type, resummation 
of a fixed order calculation includes contributions from higher order matrix elements that 
contribute at the same order in log counting. The fact that resummation can be performed 
implies that, unlike the case of clustering, there is a relationship between the logarithmic 
contributions of matrix elements at all orders. Successful resummation of non-Abelian 
NGLs would rely on this structure. Currently the leading NGLs have been resummed in 
the large N c limit, where such all orders relationships between matrix elements are used 
[42]. 

It is expected that NGLs will contribute at NLL in the exponent at all orders. For 
the anti-kT algorithm, the measurement function for the dijet mass distribution factorizes 
into single gluon constraints, meaning that only correlations in the matrix element give rise 
to the NGLs in this observable. The challenge in resumming such NGLs can be phrased 
in terms of understanding how the matrix elements that give rise to the higher order 
NGL terms are related to the lower order matrix elements. This is well understood for 
observables such as thrust, as discussed above, but no such framework exists for NGLs. If 
this framework can be developed, it may prove applicable for resummation of a wide range 
of observables. 

6. Conclusions 

Clustering due to jet algorithms can introduce correlations into soft gluon phase space con- 
straints and ruin the simple picture of resummation, spoiling even Abelian exponentiation. 
This effect has been known for some time [43] and studied in specific cases [10, 11, 12]. 
On the other hand, NGLs are traditionally thought of as arising from correlations in the 
soft gluon matrix element for measurements in separate regions of phase space, occurring 
in the non-Abelian terms first at o? s CfCa ln 2 (Q/c/A). Considering clustering logarithms 
in this light, we found it natural to interpret them as a class of NGLs. Like traditional 
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non-Abelian NGLs, they arise in the limit of gluons becoming soft rather than collinear. 
They are not associated with UV divergences, and therefore do not effect the anomalous 
dimension of the matrix elements in the factorized cross section and consequently are not 
resummed by standard RGE. 

Correlations in the measurement function arise in a wide range of jet observables. We 
considered a straightforward observable to highlight this effect, the e + e~ —> 2 jet cross 
section with measurements of the total jet mass p and the energy outside the jet A for 
the anti-kx, C/A, and bp jet algorithms. We worked in the limit {Qp, A} <C Qyfp <C Q, 
where SCET provides the relevant description, and took the jet size R S> J~p. In this 
regime the soft function contains all of the nontrivial jet algorithm dependence. In Sec. 2, 
we showed that the anti-kp jet algorithm obeys Abelian exponentiation and gave its all- 
orders form. This is because the anti-kp algorithm does not introduce correlations in soft 
gluon phase space constraints, which has been shown to be true to leading order in SCET 
power counting [36]. We also examined the violation to Abelian exponentiation due to 
the C/A and kp algorithms, using constraints from the boundary clustering phase space 
regions to show that for n final state gluons, the leading clustering logs are at least of order 
C«ln"Qp/A), which is NLL. 

In Sec. 3, we calculated the Abelian soft function at 0(a1) up to constant terms for 
the C/A and kp algorithms, including the full R dependence of the coefficients of the 
logarithms. We verified our results of the leading and next-to leading Abelian NGLs are in 
agreement with the distribution da /dp from the Monte Carlo program EVENT2. Taking 
the small R limit for the leading clustering logarithm for C/A, we reproduce the result of 
[12]. Analogous clustering NGLs exist for other factorizable jet algorithms, with calculable 
coefficients. 

We outlined the properties of clustering NGLs in Sec. 4 and compared these to tradi- 
tional NGLs, showing that their key features are the same when examined in the context 
of SCET. In Sec. 5, we discussed the all-orders structure of NGLs. We argued that at each 
order a™, there is a new correction term ASM that comes from clustering effects of n soft 
particles contributing at least at NLL. We show that this is the case in later work [22] 
by deriving an all orders form of the measurement function, which allows us to prove the 
all-orders form of the Abelian soft function. This has serious implications for resummation. 
Since each correction term AS^ for n > 2 contributes at least at NLL and the coefficients 
for different n appear to be unrelated, it seems impossible to resum the Abelian contri- 
butions to the cross section at NLL using traditional log counting in the exponent of the 
distribution. We also discussed the prospect for resummation of traditional NGLs. These 
NGLs are based on the perturbative structure of eikonal matrix elements in QCD. Since re- 
summation exploits relationships between the logarithmic structure that arises from these 
matrix elements, it is likely that a similar relationship exists for non-Abelian NGLs. 
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A. Calculation of 0(a 2 s ) Abelian Non-Global Logarithms for the C/A and 
kx Algorithms 

In this appendix we give the details of the 0{a 2 s ) calculations of the Abelian non-global 
logarithms in the C/A and k/r dijet distributions in the sum of the jet masses, Qp, and 
energy outside the jets, A. To determine the violation to Abelian exponentiation we focus 

(2) 

on in Eq. (2.13), which corresponds to the region of phase space where the algorithm 

clustering differs from the anti-kx algorithm and gives rise to a change in the cross section 
from anti-kx- 

The Abelian a 2 s soft function correction term due to clustering is given in Eq. (3.2). 
Recall Abelian contributions (in this case to a 2 Cp) with virtual gluons are scaleless in 
pure dimensional regularization and so we do not consider them here. The phase space 
constraints and matrix element only depend on the light cone coordinates ,k^ , and 
the angle between k\ and k2 in the plane transverse to the jet axes. Therefore we can 
perform integrals over the other momentum components using the on-shell conditions. In 
pure dimensional regularization with d = 4 — 2e dimensions in the MS scheme, the soft 
function correction term is 



(2),. k \ „/^C f V(MW frrgra-q, /" _ 2< 




dk+ dk^ dk+ dk^ {k+k^k+k^)- 1 - 6 AM^ip, A; k x , k 2 ) . (A.l 



o 



The calculation of the correction terms are simplest using the coordinates of the energy 
and angle of the gluons, 

kf = k°(l - cos 9i) , kj = k°(l + cos Oi) , => dkfdkr = 2 kf dk^dcos Oi , (A.2) 

where the angle Oi is measured with respect to the n direction. Since the clustering of soft 
gluons across the jet boundary is free of collinear divergences with the jet direction, by 

(2) 

IR safety of the jet algorithm, divergences will arise in the measurement function AAd^ 
in the and k® integrals. It is simpler to extract these divergences if we change the u 
coordinates to ut and z, 

k-f- = k-^ ~\~ A^2 j z = ~~q ~q , = r > dk-^dJv2 ~~ dk^ dz . (A. 3) 

A.l The Measurement Functions for the Correction Terms 

(2) 

The terms that contribute to AA^^ are shown in Fig. 1 and given schematically in 
Eqs. 3.4, 3.5, and 3.6. We now express the measurement functions in terms of the variables 
that are convenient for the calculation of the correction terms, {z, k® , cos Ox, cos 82, cos cp}. 
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For C/A, we can write the constraints @J^ A and Q®® A purely in terms of cos 9\, cos 02, 
and cos 4>: 

©C/a 01)^2,0) = 6>(cos6>i < cos 6>i 2 ) O^cosR < cos^i < cos ^ 

- COS R < COS #2 < cos R 



egJ? A (i?, 6»i,6»2,0) = 6^ -cosi? < cos 6*1,2 < cosiij 9[^cosR< cos 9 12 j , (A.4) 

where the angle #12 is given by 

cos 6*12 = sin 9\ sin 6*2 cos cj> + cos 9\ cos 9 2 , (A. 5) 

and we have restricted 9\ and 6*2, the angles of k\ and k 2 with respect to the n-jet direction, 
to the ranges where clustering can actually take place. Note that configurations related by 
the interchange of gluons 1 and 2 or the two jets are accounted for by symmetry factors in 
AXgg, as in Eq. (3.4). 

For Ict, the constraints on clustering depend on the relative energies of the two gluons. 
We note first that the out-out constraint &® 1S ^ ne same as C/A, 

© k ° T ° = ©8/°A , (A.6) 

since the constraint only requires that the two gluons are initially outside the jets, are 
separated by an angle 9\ 2 < R, and are clustered into the jets. The in-out constraint 
differs, however. If gluon 1 is in the jet and gluon 2 is out, and they cluster before gluon 1 
clusters with the jet, then the constraint is 

mm(E 2 ,E%)(l - cos 9 12 ) < E%(1 - cos 61) . (A.7) 

Clearly if E\ < E 2 , then this produces the same constraint as C/A, that 9\ 2 < 9\. If 
E 2 < Ei, the constraint on #12 is now looser: 



E 2 

1 - cos 6*i2 < -=| (1 - cos 6*i) . (Ai 
Ei 



2 

This tells us that any time the two gluons will cluster with C/A, they will cluster with 
kx, and kp has an additional region of phase space with #12 > 61 and E 2 < E\ where the 
gluons will cluster. It is useful to write the in-out constraint for kx, therefore, as 

eg = ©c°/a + ©k (A.9) 

where 

Q l g = 9(cos9 X2 < cos#i) 9(z 2 (l - cos0 12 ) < (1 - z) 2 (l-cos9 1 



cosR < cos 6*i < lj 6"^ - cos it! < cos6> 2 < cos Rj . (A.IO) 

The constraint that E 2 < E\ is a consequence of the constraints 9\ 2 > #1 and the constraint 
in Eq. (A. 8). Note that 0^ is independent of k®. 
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Soft divergences arise from each of the in-out terms -M.f at 1 /e 2 and the out-out terms 
M °\g at 1 / e - These diver 

gences come from the integrals over k®, k® (or equivalently k® , z). 
Whether the combined soft gluon, kf, is in or out of the jet is given by the angle 6t of kt 
with respect to the n-jet direction. If 9t < R, then the clustered gluons are in the jet, and 
otherwise they are out. This angle can be written as 



cos 2 8 t 



((1 — z) cos 6\ + z cos 82) 



(A.ll) 



1 - 2z(l -z){l- cos 9 12 ) ' 
Note that here 9t is independent of k®. We can now express the in-out contributions as 



M 



10 



Q 6 (cos 2 R < cos 2 6t) 6[Qp 



A: 



(1-2)(1-cos0i) + z(1-cos0 2 ) 5(A) 



Mf = 9 (cos 2 6 t < cos 2 R) 6(p) 5 (A - k° t ) 

Mf = -Q5(q P -k° t (l-z)(l- cosfli)) 5 (A - k° t z) 



and the out-out contribution as 
M°° = Q9 (cos 2 R < cos 2 9 t 



M 



p 
00 



6[Qp 



(1 - z)(l - cosfli) + z(l - cos(9 2 ) 



-9 (cos 2 R < cos 2 t ) S(p) 5 (A- k° t ) . 



(A.12) 

5(A) 
(A.13) 



Note that these parts of the measurement function are independent of the algorithm, and 
apply to both C/ A and kx • We find it is easiest to start with the calculation for C/ A and 
then tackle the kx calculation. 

A. 2 Soft Function Correction Term for C/A Dijet Mass Distribution 

We will evaluate the soft function correction for each in-out and out-out measurement 
contribution in Eq. (3.4), where ASj° c ^ A and AS®q, a correspond to A4j° and M®° 
respectively. The a 2 Abelian soft function correction term due to C/A clustering is then 
expressed as 



ASg) A (p,A) 



^^p°C/A + AS^C/A + A^i^c/A 



+ 



(A.14) 



Our general approach for each term will be to perform the k^ and z integrals to extract the 
divergences in each term and then to perform the remaining integrals numerically. We will 
find that each term can be expanded in terms of finite numeric integrals that have a common 
form. This will allow us to show analytically that the divergences in the correction terms 
cancel and yield a simple form for the result, where we need only numerically calculate 
coefficients of logs. 

A. 2.1 In-out Contribution 

We begin by considering the in-out contributions, where each term AS 1 ? can be written 



as 



7T 



a s C F \ 2 (e^) 2£ /V/ 2 r(i 



r(i - e) 




d cos <J) d cos 0i d cos 62 



(A.15) 



X (1 - cos 



: (i - cos 2 ^r 1 -^! - cos 2 e 2 )- 1 - e 8 ef /A (R, e 1 ,e 2 , <p) i\° 
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where ij contains the integral over k® and z from which we shall extract the soft diver- 
gences 



If = ^ 



dkU dz{k^ e {{l-z)zr^Mf. 



(A.16) 



'0 JO 

For c / A , the anti-kx term, this is straightforward: 



If = ^ 



dkU dz{^y^{{l-z)z)- 1 ^Mf 



Jo 
4e/n-2e „-l-2e A -l-2e 



-IjFQ-^p- 1 -**- 1 - 2 * (1 - cosi?) 



2r 



i + 2Hn|^ COS |)+e 2 c^ 
1 — cos R 



(A.17) 



(2) 

Here c a only contributes to the finite, p and A independent terms in the soft function. 

For the AS l ® c ^ A and AS^ C ^ A terms the z integral is more complex since the constraint 
on 0f n depends on z. Singularities in the z integral exist at z = and z = 1, corresponding 
to either gluon becoming soft. The constraint on z is quadratic if we write it in terms of 
cos 2 0t and cos 2 R using Eq. (A. 11). Physically, we can see that when z = 0, the gluon out 
of the jet is much softer and we will have 9% < R. When 2 = 1, the gluon in the jet is much 
softer and we will have 9t > R. This means that given cos 9%, cos 02 , and cos^> there will 
be some zq with < z$ < 1 that is the crossover point where Q±{z = Zq) = R. So we can 
rewrite the constraints on 9t as 



#(cos 2 R < cos 2 6 t ) = 0(0 < z < zo(cos #i, cos 62, cos (/>)) , 
9 (cos 2 9t < cos 2 R) = 6>(zo(cos 9\, cos 92, cos (f>) < z < l) . 



(A.18) 



We can solve for zq by solving the quadratic equation, az 2 + bz + c = 0, when cos 2 9m 
cos 2 R, taking the root 4 



z = ^{-b- Vb 2 - 4:ac) , 



(A.19) 



where 



a = (cos 9\ — cos 92) 2 — 2 cos 2 R(l — cos #12) , 

b = —2 cos 0i (cos 0i — cos 02) + 2 cos 2 i?(l — cos #12) , 

c = cos 2 0i — cos 2 R . 



(A.20) 



Performing the integral over z analytically is still challenging. Instead, we note that the 
integral over k® produces a single divergences, and therefore to obtain the divergent and 
finite terms in the z integral we only need to expand the z integrand over z to 0(e). For 
the ASyP C ; A term, corresponding to the case where the gluons are clustered with C/A and 

4 We have argued that a single root exists in the range 2 £ (0, 1). One can show that this is the smaller 
root if a > and the larger root if a < 0. In either case, this corresponds to choosing the negative sign in 
the quadratic formula. The root with the positive sign corresponds to a value of z outside of the physical 
range. 
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end up outside the jet, the integrals are 



rIO _ .M 
1 A — M 



dtf / dz 11$ 



Ox -1-46 



4e a -l-4e 



2c 



/x A 



((1 - z)*) -1-2 * #(z < z < 1) <5(p) <J(A - A;?) 



1 - z 



,2„(2) 
-A 



(A.21) 



where as with AS^ C/A , c A ^ will only contribute to the finite, p and A independent parts 



of the soft function. Similarly for AS* C / A , 



dz (jfep) 1 4e ((1 - z)^)- 1 -^ 0(0 < z < z ) ^(A) 



x Q5[Qp-k° t (1-z)(1-cos0i) + *(1-cos0 2 ) 



(A.22) 



// 4e 



2c 



Q~ 4e <5(A) p" 1 " 46 (1 - cosi?) 



l + 2e 



ln^ + 21n 1 ~ C ° Sgl )+ e M 2 ) 
1-zo 1-cosRl 



(2) 

As above, c p will only contribute to the finite, p and A independent parts of the soft 
function. 

Now we have extracted all of the divergences, and only finite integrals remain. At this 
point, we assemble the results of the in-out contribution to the soft function, and we will 
find that we can simplify the result significantly before evaluating the remaining integrals 
numerically. For each term ASj q, a expressed as in Eq. (A. 15), after doing the k® and z 



integrals we have 



jIO 
rIO 

rIO 



1 



4e a -l-4e 



// e A 



1 + + 



2c 

~Q~ ie p~ l ~ ie (\ - cos i?) 4e <5(A) [1 + e C « + e 2 C ( 2 )] 



(A.23) 



where 



-"a 



1 — COS $1 
1 — cos ' 



4 1) =21n 

1 - z 



,(1) 



-21n^+41n 1 - C ° S01 



1 - 2b 



1 — cos R 



(A.24) 



We can see that for each term, the integral over the leading order divergence is the same. 
At the next order, the c^p change the numeric integrals, but they are related since 

-2ci 1 ) + 4 1) +4 1 ) = (A.25) 



We have added R dependent terms into and (and compensating terms into the 
prefactors, which will change the arguments of the NGLs) to ensure that the remaining 
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integrals are finite as R — > 0. We will write the soft function correction in terms of the 
remaining integrals. Define 

\ -l 



A 



i,C/A 



(R) 



(e») 



2c 



TT 



r(i- e) 2 i r(i- e ) 



d cos (j) d cos 8\ d cos 62 



x (1 - cos 2 </0~ £ (l - cos 2 0i)- 1-£ (l - cos 2 6 2 )~ 1 - e c 



,(°) 



® 1 c/a( R '° 

1(0) 



Is ^2) 



(A.26) 



,(J) 



1 for all i and we just call this integral A c j A (R). Since we 



for j = 0, 1, or 2, where c, 
will find the total result is finite, we will be able to set e = in the integrals and perform 
them numerically. Using Eq. (A. 25), we have 



-2A 



(i) 

a,C/A 



+ A 



(1) 

A,C/A 



+ A 



(i) 

P.C/A 







(A.27) 



In terms of these integrals, the in-out contribution to the soft function correction term is 
very simple: 



1 



7T / Q(l — cos R) 



x < 2Cg) A (R) 



P 



1 — cos R 



P 



(i) 

IO.C/A 



(2) 



1 — cos i? 



+ 5 
A 



p 



1 — cos R 



6 



Cr 



P 



P 



VI -cos 12 

where jCqOe) an d C\{x) are distribution functions. 



1 — cos R 
A 

Q 



p 

1 — cos 1? 



Co 



A 



C (x) 

The NGL coefficients are 

-.(0) / 
y C/A l 



9(x) 
x 



C^x) 



9{x 



lnx 



C 



c$] a (r) = 1 -a<-°Hr), 
1 

2 



(A.28) 



(A.29) 



(A.30) 



(i) 



IO.C/A 



(1?) 



.4 



C 1 ) cm _ 4W 

ci.C/aW p,C/A 



(R)) , 



We will not explicitly calculate C^c/A 

(12). 

We can see that the divergences in AS 10 have cancelled, which is a consequence of 
Eq. (A.27). This means we can set e = in the ^^/A integrals. Note that we have shifted 
the argument of the single NGL to include a factor of 1 — cos R to give the same argument 
as the double NGL. This change simply alters the constant p and A independent terms. 
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A. 2. 2 Out-out Contribution 



Next we consider the out-out contribution to the soft function correction term, represented 
in Fig. 1(c). Each term AS®° can be written as 

-l 



A cOO 



asC F y ( e ^r h^ni-e) 

r ( i- e) 2 I r(i- e ) 



d cos (i cos 6\ d cos #2 



(A.31) 



X (1 - COS 2 0)~5" e (l - COS 2 - COS 2 02)^^4 G%? A (R, #1, 2 , 



OO 



where as before I® contains the integrals over k® and z from which we shall extract the 
soft divergences: 



l 



OO 



-l-2e 



OO 



(A.32) 



Since both gluons are outside of the cone of radius R, as z — > or 1 the combined pair k t 
will be outside the jet. Therefore, there are no divergences in the z integral, and the only 
divergences in the out-out contribution arise when both gluons are soft (k° = 0). As in the 
in-out case, the constraint cos 2 R < cos 2 8t can be expressed as a quadratic function of z, 

9(cos 2 R < cos 2 9t) = #(zi(cos #1, cos #2, cos </>) < z < zq(cos 9\, cos 62, cos <f>)) , (A. 33) 

where 

zi = ^-(-b+Vb 2 -4ac) , zo = —(-b- \A 2 -4ac), (A.34) 

_ ( I —(I 

and a,b,c are defined in Eq. (A. 20). Since zq and z\ are independent of k\ we can carry out 
this integral and extract the divergent contributions. The contribution from C/A clustering 
for this configuration is 



dtf 



dz (k 



\-l-4e 



((l-z)z) 



-l-2e 



x5[k° t 



Qp 



(1 - z)(l - cos6»i) + z(l - cos( 



<5(A) 6>(zi < z < z ) 
\ ((1 - z)(l - cos6»i) + z(l - cos(9 2 )) 



■k 



1 //i 



k 



-l-4e 



In 



(1 - zt)z 



+ ec. 



(2) 

00,p 



(A.35) 



(1 - zq)zx i 

where we can set the e = in the integral over z, since we have argued it is finite. The 
corresponding terms from anti-kx for the out-out configuration give 



I 



OO 
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dk° f 1 


Jo 
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dz (fc°) 1 4e ((1 - z)z)~ l ~' ze 5(p) 0{zi <z<z ) 5(A - jfe°) 



-l-2e 
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-l-4e 
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(1 - 2q)zi 



+ ec 
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OQ,a 



(A.36) 



We define the remaining angular integrals in the calculation of the out-out soft function 
to be 



48>.i(*) 



rfi 




<i cos d cos 0i d cos 02 (1 — cos 2 



x (1 - cos 2 fli)" 1 -^!- cos 2 tfs)- 1 -^ e°° A (i?A, 02,^)^, , 



(A.37) 
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for j = 1, 2 and i = p,a where 



In 



(1 - Z])zq 



-oo,i "^(1-^1. 
The out-out soft function correction term then takes the simple form, 



(A.38) 



a cOO _ a oOO , a oIO 
^C/A - ^p.C/A + ^>a,C/A 
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+ C, 
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00,C/A 
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4 A 

cosRj \Q 



(A.39) 



1-cosRj °\Q 



The coefficients are 



a 



(i) 

OO.C/A 



C oo,c/a(*) = ^oo(^) " «»*) + 4& >p (/l) - 4& >a (fl) 



(2) 



l (2) 



(A.40) 



As with the in-out contributions, since the divergences have cancelled we can set e = and 
evaluate the remaining integrals numerically. We do this in the following subsection. Note 
that we have shifted the argument of the out-out NGL to equal the argument of the in-out 
NGLs, which affects only the constant p and A independent finite terms. 

A.2.3 Total Soft Function Correction Term for C/A 

Combining the results from the previous subsections, we find the soft function correction 
term from clustering is 



\ c( 2 ) ( \\ I a sCF 
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Cr, 
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Cr 



+ ( C 10,C/a( R ) + ^00, C/A 



2sin 2 #) £ °(Q 



+ ( c io,c/a( r ) + c io,c/a( r )) 5 
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2 sin 



2 R 



(A.41) 



Converting to the singly differential distribution in p and the cumulant as follows, 

S{p) = dA S(p, A) S S = r dp S(p) , (A.42) 

Jo Jo 
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the p-dependent terms in the correction to the singly differential distribution and cumulant 
are 

2 



a s CF 



7T 



2C c/a( r ) 



+ c 



,(1) 



(1) 



IO,C/AV x V 1 ^00,C/A 
(2) , m . ~(2) 



£o(p)-^(p)ln 



2A C . 2 R 
IT" 11 2 



+ riO.C/A^) + ^IO.C/A 



(A.43) 



A 4 2 C/A- (^) {c$)aW ln ' 2A Q S ^ | + ( C IO,C/A(^) + 4o,C/aW) 111 2A Q ^2 R 

+ constant I , (A. 44) 

where the constant are two loop finite terms that we do not determine. Since the total 
result is finite, we can set e = in the A^°\A^ and ^4oo integrals and perform them 
numerically. These integrals reduce to 

A iC/A ( R ) = ~ T d cos <t> d cos i d cos ^ T , (1 I fl Tn^ 2 27m 0I c°/a ^i.«2, <A) c?' } 

'• U ' A 7T 7 — 1 (1 — COS — COS O2) ' 

4. f 1 (1 — COS 2 0) ^ 

^oo(-R) = - / d cos calicos 6*i d cos 8 2 T , 9 „ W1 5^-r®c/A ^2, 0) coo ■ 

7T y_ 1 (1 — COS^ <7l)(l — COS^ #2) 7 

(A.45) 

These determine the coefficients of the clustering NGLs, which we plot in the main text in 
Fig. 7. 



B. Soft Function Correction Term for the k T Dijet Mass Distribution 

The soft function correction term for kx also comes in two parts, an in-out contribution 
and an out-out contribution. As discussed previously, the out-out contribution is the same 
as C/A, and was calculated above. The in-out contribution for bp differs, and can be 
written as the in-out contribution from C/A plus an additional term, as in Eq. (A. 9). The 
clustering constraints for this additional term, 0^*, are defined in Eq. (A. 10). We will 
calculation the soft function correction term, AS^*, that corresponds to 0^*. Just as for 
C/A, this has the form 

AS™ = A5i° A + ASl° A + AS*° A , (B.l) 

with 

ASl° A = (^) 2 f0^ ( ![ ^~ir ) //W^os0 lC Zcos0 2 (B.2) 
x (1 - cos 2 (p)-^~ e {l - cos 2 tfi)- 1 "^! - cos 2 62)^8 e^ SA (R, 01, 2 , <f>) 4, A , 
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where 6^ A is equal to Q 1 ^ with the z-dependent kx clustering constraint removed, 
©PS,A = cos ^12 < cos0i) 9^cosR < cos 0i < l) - cos -R < cos 6*2 < cos-R) , (B.3) 

and lj° A contains the integrals over and z from which we shall extract the soft diver- 
gences, 

4°a = ^ r < c dz (k*)- 1 ^ ((i - z)z)-™ 

Jo Jo 

x 0(z 2 (l - cos 0i 2 ) < (1 - z) 2 (l - cosfli)) . (B.4) 

The clustering constraint that we extracted from Q l A and put into lj° A can be written as 
a constraint on z: 

e[z 2 {\ - cos0i 2 ) < (1 - z) 2 (l - cosfli)) = 9{z < z kT ) , (B.5) 

where 

(^^Hf)" 1 - 

calculation of terms is the same as the 7? terms for C/A with 
r . This is a simple change to these calculations, and the results 
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where 
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4e A -l-4es 



(l-z kT )z 5 ' A 



0(z o < z kT ) 



(B.8) 
(B.9) 

(B.10) 



Note that 7*°^ has p and A dependence that is not simply a distribution in each variable, 
since there is a constraint on the relative values. This means it is difficult to expand 
the contributions out in terms of distribution (plus) functions. Therefore we will only 
determine the double cumulant soft function for the kx algorithm, which is much simpler. 
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We are free to move to the cumulant before performing the remaining finite integrals, and 
for each contribution we obtain 



r io,s 

a, A 




sin ^ sin e -f Qp c 




sin 2 f 



sin 2 f 



Pc 

1 — cos R 

> 



< 



2A C sin 2 f 



-2e 



~ 2t ( sin \ sin ' 
{ sin 2 f 



sin % sin Qp c 



2A C sin 2 f / V 1 - cos R 



x 



l + 2 e lnf i- C ° S M +^i 2) A 
\l-cosRJ a - A 



;) ^(l-cosi?)- 2 ^ 



(B.ll) 



rIO.E 



4c 



1 

8? 



Pc 



1 — COS i? 

l + 2e| -In 



-4e 



+ 2 In 



1 — COS #1 

1 — cos R 



+ e 2 c (2) 
+ e C P ,A 



(B.12) 



'A, A 



/i 



k 



1 /A, 



-4c 



If we add the cumulants, then we find 
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(B.14) 



Note that all the divergences cancel. The constant depends only on R, and we do not 
calculate it. Therefore, as with C/A, we can set e = and define the coefficients in terms 
of finite integrals which we will evaluate numerically. 
The coefficients are related to the integrals 
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d cos 4>d cos 0!d cos 9 2 {1 - cos 2 0) 2(1 _ cos 2 6>i) 1 (l-cos 2 2 ) 1 



x ^ G^s A (R, Oi, 8 2 , 0) c%> , 



Si) 



(B.15) 
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where 



mi 1 / sin -if sin -s 4 , , , 

' " " J ' 2R < \* I . (B.16) 



sm ¥ 

012 



In 




.(1) _ g / Shl 2 Sln "2^ ^ QPc 1 l f Z k T 1 ~ COS 



1 — Zk T 1 — cos R 



' zp(l - z kT ) 
^k T (l - z ) 

Then the soft function cumulant AS^ A (which is the cumulant of S^P defined in Eq. (B.l)) 



is 

f J \.2A c sm ^ 



AS^ = A%> In 2 ^ c „ B + A% In a R + constant. (B.17) 
A \2A c sm 2 § J A ^Acsi^fy 

Finally, we can add these terms to the soft function cumulant for C/A and obtain the result 



for k 



T- 



A41 = (^) {<>(«) ^ ^ + (cg lT ,H, + eg- tT («>) , ^ 

+ constant I , (B.18) 



where 



rr(o) _ r (0) 4 (0) 
°k T - °C/A A ' 



n m _ r rW ,4(1) 
°IO, k T ~~ °IO, C/A A J 

C oo, k T = C'oo, C/A • (B.19) 

Note that the in-out k^ coefficients implicitly depend on the values of p c and A c through 
A^' 1 ^ 1 . In particular, they depend on the value of the argument of the NGL, Qp c /2A C sin 2 ^ . 
From Eq. (B.16), we can see that when this quantity is larger than 1, the p c and A c 
dependent constraint will always be satisfied, since 6\ < R and 812 < R- As the NGL 
argument decreases, it will shrink the phase space over which the coefficients contribute, 
reducing their value. These coefficients are plotted as a function of R in Fig. 7 for the 
range Qp c /2A c sin 2 f > 1. 
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